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ABSTRACT 

Bow shocks are ubiquitous astrophysical phenomena resulting from the supersonic passage of an 
object through a gas. Recently, pre-transit absorption in UV metal transitions of the hot Jupiter 
exoplanets HD 189733b and WASP12-b have been interpreted as being caused by material compressed 
in a planetary bow shock. Here we present a robust detection of a time-resolved pre-transit, as well 
as in-transit, absorption signature around the hot Jupiter exoplanet HD 189733b using high spectral 
resolution observations of several hydrogen Balmer lines. The line shape of the pre-transit feature 
and the shape of the time series absorption provide the strongest constraints on the morphology and 
physical characteristics of extended structures around an exoplanet. The in-transit measurements 
confirm the previous exospheric Ha detection although the absorption depth measured here is ^50% 
lower. The pre-transit absorption feature occurs 125 minutes before the predicted optical transit, a 
projected linear distance from the planet to the stellar disk of 7.2 Rp. The absorption strength observed 
in the Balmer lines indicates an optically thick, but physically small, geometry. We model this signal as 
the early ingress of a planetary bow shock. If the bow shock is mediated by a planetary magnetosphere, 
the large standoff distance derived from the model suggests a large planetary magnetic field strength 
of Beq=28 G. Better knowledge of exoplanet magnetic field strengths is crucial to understanding the 
role these fields play in planetary evolution and the potential development of life on planets in the 
habitable zone. 


1. INTRODUCTION 

Hot Jupiters (HJs) are Jupiter-mass planets that are 
orbiting within ^10 stellar radii of their host stars and 
have orbital periods of a few days. The large insolation 
of these planets results in extreme heating of their bound 
atmospheres and inflated radii (Burrows et al. 2007). 
These objects provide a view into planetary atmospheric 
physics that do not occur in our solar system. Some HJs 
are believed to be losing mass via hydrodynamic blowoff 
(Vidal-Madjar et al. 2003). The close proximity of a HJ 
to its parent star can result in star-planet-interactions 
(SPIs) (Cuntz et al. 2000; Shkolnik et al. 2005) and, po¬ 
tentially, the formation of bow shocks around the planet 
as it moves supersonically through the stellar wind or 
corona (Vidotto et al. 2010). If the bow shock around a 
HJ is mediated by the planetary magnetosphere, absorp¬ 
tion by material compressed by the shock can be used to 
estimate the strength of the planetary magnetic field. 
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One of the closest of these HJs to our solar system, 
HD 189733b (Bouchy et al. 2005), has been studied ex¬ 
tensively due to the relative brightness of its host star. 
Various molecules and atoms in both the bound and un¬ 
bound atmosphere of HD 189733b, including (but not 
limited to) GH 4 (Swain et al. 2008), H 2 O (Grillmair et 
al. 2008), Na I (Redfield et al. 2008), O I (Ben-Jaffel 
& Ballester 2013), and neutral hydrogen (Jensen et al. 
2012; Lecavelier des Etangs et al. 2010, 2012), have been 
detected using transmission spectroscopy, where stellar 
photons are absorbed by the planetary atmosphere while 
the planet transits its host star. 

In addition to in-transit absorption, hints of pre-transit 
absorption have been observed in HD 189733b (Ben- 
Jaffel & Ballester 2013; Bourrier et al. 2013) and WASP- 
12b (Fossati et al. 2010) and provided early indications 
that material exterior to an exoplanet’s atmosphere could 
have sufficient optical depth to cause an absorption sig¬ 
nature. The WASP-12b result, a single 2-sigma absorp¬ 
tion of a 41 A region in the near-UV, was interpreted and 
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modeled as a transiting bow shock (Vidotto et al. 2010; 
Llama et al. 2011). Bourrier et al. (2013) measured pre¬ 
transit absorption in Si III 1206.5 A for one of two ana¬ 
lyzed HD 189733b transits. Ben-Jaffel & Ballester (2013) 
report two pre-transit UV absorption measurements for 
HD 189733b, observed in C H, but not O I, and they 
model the result as a transiting bow shock (Ben-Jaffel & 
Ballester 2013; Llama et al. 2013). 

Here, we present a time-resolved pre-transit signal de¬ 
tected in the high-resolution transmission spectra of the 
Balmer lines Ha, H/3, and Hy of HD 189733b. The ob¬ 
servations and data reduction procedures are given in 
section 2. The transmission spectrum analysis and our 
procedure for estimating uncertainties in the absorption 
signal is presented in section 3. In section 4 and section 5 
we discuss the absorption signals and the effect of star 
spots and active regions. The bow shock and exosphere 
models are described in section 6. We discuss the plan¬ 
etary magnetic field strength estimate in section 7 and 
provide a brief conclusion in section 8. 

2. OBSERVATIONS AND DATA REDUCTION 

Two transits of HD 189733b were obtained on the two 
half nights of June 3 and July 4, 2013 using HiRES on 
Keck I (Vogt et al. 1994). The approximate resolving 
power of the observations is i?~68,000, or 4.4 km s“^ 
at Ha. The B2 decker was employed which has a slit 
size of 7.0" X 0.57". Exposure times ranged from 3 to 
5 minutes. The average signal-to-noise of the extracted 
spectra is 400 at Ha, 180 at H/3, and 120 at Hy. 

The observations from June 3 begin at first contact and 
so add no significance to the pre-transit detection. In 
addition, the cadence of the June observations was twice 
as long, due to intermittent terrestrial clouds and haze, 
providing less detail on the shape of the transit. Eor these 
reasons we do not include the June data in the analysis. 
We note that the in-transit absorption values measured 
for the June data are similar to the values measured from 
the July data. 

The data were reduced using the HiRES Redux pack¬ 
age written by Jason X. Prochaska^. Standard reduction 
steps were taken including bias subtraction, flat fielding, 
and the removal of cosmic rays and hot pixels. The spec¬ 
tra were extracted using a 6.6" boxcar. All images were 
examined manually for order overlap in the blue chip. 
HiRES Redux also performs 2D wavelength solutions us¬ 
ing Th-Ar lamp exposures taken at the beginning of each 
night. The residuals for the wavelength fits are ^0.05 
pixels, or ~0.02 A, in all orders. When applied to an 
individual observation, all wavelength solutions are cor¬ 
rected for Earth’s heliocentric velocity and for the radial 
velocity of HD 189733. This places the spectrum in the 
rest frame of the star. 

Telluric absorption is removed from the Ha spectra 
using the telluric fitting program Molecfit (Kausch et al. 
2014). We first remove the strong Ha line absorption 
from a telluric standard. The telluric absorption in the 
standard is then fit using Molecfit. The best telluric fit 
is then shifted appropriately in wavelength, scaled, and 
subtracted from each individual observation. This pro¬ 
cess results in transmission spectra with residuals in the 
telluric lines of ~0.25%. Telluric absorption from Mauna 

^ http://www.ucolick.org/~xavier/HIRedux/ 


Kea is negligible in the H/3 and Hy orders. Telluric re¬ 
moval is not performed for these lines. 

3. TRANSMISSION SPECTRUM AND EMC ANALYSIS 
The transmission spectrum is defined as 


where E/ is a single observation and Font is the master 
post-transit spectrum. In order to obtain the final trans¬ 
mission spectra all observations must be normalized and 
aligned in wavelength space. Changes in spectral resolu¬ 
tion due to temperature or instrument flexure are mini¬ 
mal and do not affect the spectra over the course of one 
half night. All spectra are initially normalized using a 
4th order polynomial. A single out-of-transit spectrum 
is then chosen as the comparison spectrum in order to 
apply a wavelength correction to all spectra. A number 
of stellar absorption lines in each spectrum are cross cor¬ 
related with the same lines in the comparison spectrum. 
We tested three different types of fits to the resulting 
wavelength shifts: constant, linear, and spline fits. The 
spline and linear fits perform no better, and in some cases 
worse, than a constant offset in wavelength. Thus the 
constant offsets are used to rectify the wavelength vectors 
for each spectrum relative to the comparison spectrum. 
This results in precisely aligned spectra where the largest 
residuals occur across 1-2 pixels at the cores of narrow 
spectral lines. Cores of identified stellar absorption lines 
in the normalized spectra are ignored when calculating 
the absorbed flux in each line so these residuals do not 
contribute to the transmission signal. Once spectra have 
been normalized and wavelength shifted they are divided 
by the master out-of-transit spectrum, a co-addition of 
the selected out-of-transit spectra. The resulting spec¬ 
trum is renormalized to remove slowly varying contin¬ 
uum shapes leftover from the division. This occurs be¬ 
cause the original normalization is slightly different for 
each spectrum. The line of interest is ignored when de¬ 
termining the renormalization so as to avoid introducing 
offsets from the true continuum. 

Eigure I shows the average transmission spectrum for 
the pre-transit (first column), in-transit (second column), 
and post-transit (third column) observations for Ha (top 
row), H/3 (middle row), and Hy (bottom row). The 
model line profiles are over plotted with solid lines (see 
section 4). The number of spectra used to create each 
spectrum is give in the bottom-right of the top-row win¬ 
dows. The fourth and fifth columns shows the results 
of an empirical Monte Carlo (EMC) process that we uti¬ 
lize to produce estimates for the absorption uncertainties 
(Redfield et al. 2008; Jensen et al. 2011; Wyttenbach et 
al. 2015). The EMC distributions are given in terms of 
Wa, the equivalent width of absorption across the line. 
W\ is calculated by integrating the line profile from —200 
to -1-200 km s“^. Values of Wx do not change by more 
than 25%, for any measurement, if the integration width 
is extended to ±500 km s“^. 

The in-transit EMC distributions are generated by 
randomly selecting 5000 subsets of in-transit spectra and 
comparing them to either the master out-of-transit spec¬ 
trum (the In-Out method) or the master in-transit spec¬ 
trum (the In-In method). The standard deviations of 
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Figure 1. The averaged transmission spectra for the pre-, in-, and post-transit spectra for Ha (top row), H/3 (middle row), and H 7 
(bottom row). The cores of identified stellar lines other than the lines of interest are masked in gray. These points are not included in 
the calculation of the absorbed flux. The colored lines are the average model line profiles at the times of the average measured absorption 
(see section 6 ). The fourth frame in each row shows the in-transit empirical Monte Carlo (EMC) W}, distributions for each line. Each 
absorption measurement is detected above the Scr level. The fifth column shows the EMC for the average pre—transit signal. 


the resulting distributions are taken to be the Icr un¬ 
certainties for the reported in-transit absorption. Each 
line is clearly detected at >3a. We note that propagated 
Poisson errors result in very similar estimates (^20% less 
than the EMC estimates) for the absorption uncertainties 
in each line. We choose to use the uncertainty derived 
from the EMC due to its ability to highlight systematic 
errors, which can dominate the measurement if not taken 
into account. 

The small number of pre-transit observations (N=4) 
does not allow sufficient combinations of spectra to be 
generated for the EMC process described above. In¬ 
stead, we compare the average pre-transit spectrum to 
all possible combinations of the nine post-transit spectra 
(iV=511). The standard deviation of the resulting distri¬ 
bution, shown in the hfth column of Figure 1, is taken as 
the 1-CT uncertainty for the measured absorption. This 
process is also used to estimate an average uncertainty 
for the individual points in Figure 3. 

As an additional test of the EMC procedure and to 
verify that we are actually measuring absorption in the 
Balmer lines, we perform the same in-transit EMC anal¬ 
ysis for three control lines of Ca I. One line is chosen each 
for Ha, H/?, and H 7 from the same spectral order that 
the Balmer line is located in. Ca I is expected to con¬ 
dense out of the atmospheres of hot Jupiter exoplanets 
and thus should not exist in atomic form in large quan¬ 
tities (Fodders 1999). If systematics from the reduction 
procedure are resulting in absorption in the Balmer lines 
we are likely to see similar artifacts in the Ca I control 
lines. The results of the control line transmission and 


EMC analyses are shown in Figure 2. As expected, these 
lines show no evidence of absorption and the resulting 
uncertainties are similar to the uncertainties measured 
for the Balmer lines. 

4. IN-TRANSIT ABSORPTION 

The second column of Figure 1 shows the average in¬ 
transit absorption prohles for the July transit. The in¬ 
transit EMC is shown in the fourth column. The ab¬ 
sorbed fluxes measured in the line prohles of Ha, H/3, and 
H 7 are all detected at >3a. The H,d and H 7 detections 
are the hrst reported measurements for an exoplanet at¬ 
mosphere. The absorption that we measure is ^50% 
lower than that from Jensen et al. (2012). This may 
be due to variability in the surface features of HD189733 
(e.g., spots and active regions). We model the in-transit 
prohles as arising in the planetary exosphere (see sub¬ 
section 6.2). The possible contribution of star spots and 
active regions to the absorption prohle is discussed in 
subsection 5.1. 

5. PRE-TRANSIT ABSORPTION 

Absorption values, which are the negative of W\, deter¬ 
mined from individual exposures are shown as a function 
of time in Figure 3. The gap in the data from ^ —115 
minutes until —70 minutes was, unfortunately, used to 
obtain telluric standards. The strong pre-transit signal 
was not anticipated. The absorption is calculated by 
integrating the hux from —200 to -1-200 km s“^. The er¬ 
ror bars represent the average standard deviation of the 
EMC process using all combinations of the post-transit 
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Figure 2. Master line profiles for the Ca I control lines. The format is the same as Figure 1. Each line is consistent with no absorption 
indicating that systematics do not contribute to the absorption measured in the Balmer lines. 


spectra. The average pre-transit signal shown in Fig¬ 
ure 1 is detected at > Str for Ha and H/3. The H 7 signal 
is detected slightly below the 3a level. 

The pre-transit absorption exhibits a few important 
properties. First, the ratio of the absorbed flux between 
Ha and H/3 is ^ 1.6 — 2.0, which is much lower than in 
the optically thin limit. This indicates that the absorbing 
material is optically thick. Indeed, the line morphologies 
of the pre-transit observations appear saturated, i.e., flat 
bottomed and broad. Secondly, the absorption abruptly 
begins to decrease between t — tmid = ~70 and —66 min¬ 
utes. This requires that the material causing the pre¬ 
transit absorption must begin to exit the stellar disk be¬ 
fore the time of first optical contact (labeled tj in Fig¬ 
ure 3). This suggests that the distribution of absorbing 
material is asymmetric with respect to the planet, a point 
that is supported by the lack of any significant absorption 
seen at similar times post-transit (in contrast to GJ 436, 
for example; Kulow et al. 2014). Finally, the integrated 
absorption never exceeds ^1.3% which requires that the 
optically thick absorbing material cover a similar fraction 
of the stellar surface. 

5.1. The Effects of Star Spots, Active Regions, and 
Stellar Activity 

Star spot crossings are not unusual for HD 189733b 
(Pont et al. 2013). The transit of a star spot or active 
region can change the relative amount of stellar light be¬ 
ing absorbed by the occulting material. The transit of 
a spot-free chord can also change the relative amount 
of stellar light being absorbed, even for specific spectral 
features, if there is a non-negligible fraction of the stel¬ 


lar surface covered by spots (Berta et al. 2011). Stellar 
activity can also mimic transit features when compar¬ 
ing spectra from different times: if the star brightens or 
dims relative to the time of comparison, this signal can 
appear as the absorption of stellar flux when in truth it 
is merely variation in the amount emitted by the star. 
For our case, stellar variability due to time-dependent 
Balmer line emission from active regions can mimic the 
spectral line absorption: if the star is more active during 
the times when the out-of-transit spectra are measured, 
the filling in of the stellar Balmer line cores will appear 
as an absorption signature for any spectrum measured at 
a time when the star is less bright. 

There is ample evidence that the absorption signatures 
we measure are not caused by stellar variability. Both 
pre- and in-transit spectra are compared to the same 
out-of-transit spectra. This requires that any absorp¬ 
tion signal caused purely by stellar variability have the 
same absorption strength in both the pre- and in-transit 
spectra. This is clearly not the case. In fact, the pre¬ 
transit measurements are significantly stronger than the 
in-transit measurements despite the fact that the exo¬ 
sphere, by way of the occulting planet, transits an effec¬ 
tively smaller stellar disk than the bow shock. If the pre¬ 
transit absorption is entirely due to stellar variability, the 
in-transit absorption would necessarily be stronger than 
the pre-transit absorption. Since this relationship is not 
seen in the data, the absorption cannot be solely due 
to stellar Balmer line variability. This means that the 
pre-transit signal must be caused by absorbing material 
occulting the stellar disk. 

An indicator of stellar activity level is the S'-index mea- 
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sured from the Ca II H and K lines (Duncan et al. 1991). 
The index, shown in Figure 4 for the July transit, has 
similar values pre-transit compared to post-transit, al¬ 
though there appears to be a decrease about halfway 
through the transit. The similar S-index values both be¬ 
fore and after transit also indicate that the pre-transit 
absorption is not caused by variable stellar emission. 

The strength of spectral features can be affected by 
the relative contributions of the stellar photosphere and 
spots or active regions (Berta et al. 2011). While the 
transiting bow shock, which is transparent to continuum 
photons, is not subject to this affect, the opaque tran¬ 
siting planet certainly is. If the planet transits a spot- 
free chord, the stellar spectrum will be weighted towards 
the spot spectrum and the relative strength of spectral 
features between in and out-of-transit observations can 
change, even at the 1 — 2% level. As a result, this effect is 
capable of producing spurious in-transit features. We do 
not believe this to be the case for our in-transit absorp¬ 
tion measurements. In general, spotted active regions 
will £11 in the Balmer line cores with emission, causing 
the line strength to effectively decrease. When the planet 
is transiting, this produces weaker absorption lines since 
the observed spectrum is weighted towards the active 
spotted regions. When compared to out-of-transit mea¬ 
surements, this will produce in-transit emission features. 
The fact that we see the opposite, i.e., strong in-transit 
absorption features, indicates that we are measuring ab¬ 
sorption through the atmosphere of the planet. How¬ 
ever, the absolute strength of the measured absorption 
is almost certainly affected. The contribution of spotted 
active regions most likely causes us to underestimate the 
strength of the in-transit exospheric absorption. This 
may be the cause of the difference between the in-transit 
absorption strength measured here and that of (Jensen 
et al. 2012), who measure the in-transit line depth to be 
two times deeper. 

The arguments given above apply to the general fea¬ 
tures of the absorption time series. It is difficult, how¬ 
ever, to determine how individual absorption measure¬ 
ments are affected by a nonuniform stellar disk. Some 
absorption measurements may well be affected by spots 
or bright regions on the stellar surface. For example, 
in between 3rd and 4th contact the absorption varies 
between the expected values from the exiting exosphere 
and much stronger absorption not accounted for by the 
model. These points could be due to the transitory occul- 
tation by the exosphere of bright chromospheric regions 
on the limb of the star. Modeling of these variable fea¬ 
tures is beyond the scope of this study but they provide 
a plausible explanation for the departure of individual 
points from the model estimates. 

6 . MODELING THE ABSORPTION 

We model the pre-transit absorption as arising in a 
thin bow shock orbiting ahead HD 189733 b. The in¬ 
transit absorption is modeled as arising in the exosphere 
of the planet. We neglect limb darkening in our model 
calculations since we are measuring the absorption across 
a very narrow (^I.O A) portion of the stellar line. This 
portion of the Balmer lines forms in a very thin layer in 
the stellar atmosphere with a small temperature gradient 
and thus suffers negligibly from limb darkening (Mihalas 
1978). 


6 .1. Bow shock model 

We construct a simple model of the pre-transit ab¬ 
sorption as arising in a bow shock orbiting ahead of the 
planet. While a bow shock is not the only possible cause 
of the absorption, it is a plausible one and has been suc¬ 
cessful at explaining the pre-transit UV observations of 
WASP12-b and HD 189733b (Llama et al. 2011; Ben- 
Jaffel & Ballester 2013). We do not attempt to model 
the full physical conditions of the stellar wind/corona 
and the resulting bow shock conditions. Instead, we as¬ 
sume a bow shock geometry and require that the n = 2 
hydrogen density assigned to the bow shock produce the 
necessary absorption. We assume that the standoff dis¬ 
tance of the magnetosphere, is equal to the distance 
of the nose of the shock from the planet. 

The bow shock model is constructed using the shock 
geometry of Wilkin (1996). A density is assigned to the 
nose of the shock and then scaled along the bow accord¬ 
ing to the scaling law 

pW = Po(^) (2) 

where r is the linear distance from the planet to the bow 
shock and a is the exponent describing how quickly the 
density changes with distance from the density at the 
nose, po- The distance r is always greater than so the 
density decreases monotonically from the nose. A 3-D 
density grid is then filled according to the prescription 
of the density law with the chosen value of a. The bow 
shock is rotated around the z-axis with angle Qsh from 
the tangent to the planet’s orbital motion. The column 
density of the material covering the stellar disk is then 
calculated. The cell size in the grid is 0.01 i?p, the same 
size as the shock thickness. We find that a cell size of 
this magnitude is required in order to reproduce the rel¬ 
atively small absorption values from the data. Larger 
cells, which necessarily result in a larger shock width, 
result in too much of the stellar disk being covered by 
optically thick material. We find values of Wx to within 
^5% of the reported values if we decrease the grid cell 
size to 0.005 Rp, indicating that the shock thickness of 
0.01 Rp is not an artifact of the resolution. 

The optical depth in the Balmer lines is approximated 
by a Doppler-broadened delta function, = Toe~^'"^^'> , 
where b = and cr„ is the dispersion of a 1-D Gaus¬ 
sian velocity distribution (Draine 2011). The optical 
depth at line center is 

fiuNiXiu 

To — ----• (3) 

iriecb 

where fiu is the oscillator strength of the transition, 
is the central wavelength, is the electron mass, e is 
the electron charge, c is the speed of light in vacuum, 
and Ni is the column density of the lower energy level of 
the transition. We use the optical depth as a function of 
velocity from line center to approximate an absorption 
line profile for each grid point using the column density 
calculated from the 3D bow shock. The final line pro¬ 
file for a single model iteration (i.e., a single time of the 
bow shock transit) is the sum of each individual line pro¬ 
file for each grid point on the stellar disk. Examples of 
the pre-transit model fine profiles are overplotted on the 
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Figure 3. Absorption as a function of time from transit midpoint for a single transit of HD 189733b. First through fourth contact are 
marked with the vertical green dashed-dotted lines. A sharp decrease in the absorbed flux can be seen from —70 min< t < —40 min. 
The gap in the data from —115 minutes to —70 minutes was used for telluric standard observations. The absorption does not appear 
post-transit. The model is shown with solid lines. The uncertainties in the absorbed flux for each individual spectrum are shown in the 
upper-left. Each bar shows the average of the standard deviations of Wx obtained from the EMC for each point. 



1 ■ 1 ■ ■ 1 ■ ■ 


- 


1 + 

+ 


I 

: + + 

+ \ 

! + 

- 


! + J 



+ ! + + 

+ 


+ + 

+ 1 

+ 


+! + 

. 

+ 

+ 

1 

L 

! 


- 


i + 

- 


j + 7 



,1 ■ 


-100 -50 0 50 100 

t-t (minutes) 


Figure 4. Ca II H and K 5-index measurements taken simultane¬ 
ously with the Balmer lines. Propagated Poisson uncertainties are 
smaller than the plot symbols. Transit contact times are marked 
with the vertical green lines. There is no direct correlation with 
the transit, although there appears to be a decrease in the activity 
index at t — tmid = 0 minutes. 

measured spectra in Figure 1 We find that 6 = 4.1 km s“^ 
is able to adequately reproduce the line profile shapes. 
The final model W\ values are calculated identically to 


the measured values of W\. We point out that the model 
seems to poorly reproduce the pre-transit Hy line profile 
and absorption values. We believe this is due to the 
blending of Hy with a nearby Fe I line which causes the 
red wavelength side of the Hy profile to show weaker ab¬ 
sorption. This can be seen in Figure 1 where the model 
line profile matches the blue side of the line very well. A 
similar effect is seen for the in-transit Hy profiles. Based 
on the line profile shapes of Ha and H/3 we have no rea¬ 
son to believe the shallow red-ward portion of the Hy 
line is representative of anything physical. 

The model absorption values are shown as solid lines 
in Figure 3 and the model line profiles are shown as solid 
lines in Figure 1. A to-scale snapshot of the system 
is shown in Figure 5. The large standoff distance and 
the observed absorption time series are representative of 
the double transit suggested by Vidotto et al. (20ff). 
This model was calculated using values of rj„=12.75 i?p, 
9sh = 15°, the angle between the planet’s trajectory and 
the nose of the shock, and po = 9x10“^° g cm“^, the 
mass density of excited hydrogen at the nose of the shock. 
A scaling factor of a = 400 is required to sufficiently con¬ 
centrate the density at the nose which is necessitated by 
the loss of absorption near t/. The model reproduces the 
main features of the observed pre-transit time series: 1. 
the approximate ratio of Ha to H/3 absorption; 2. simi¬ 
lar absorption values at —125 and —70 minutes; and 3. 
the steep decrease in the absorption immediately after 
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t — tmid = —70 minutes. The model also predicts an ab¬ 
sorption maximum at t — tmid= —85 minutes and ingress 
of the bow shock at t — tmid= —175 minutes. 
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Figure 5. To-scale projections of the planet and bow shock in 
the orbital plane (top panel) and the view from Earth (bottom 
panel) for time t — t^id = — 50 minutes, soon after the bow shock 
has exited the stellar disk. The planet is in between first and 
second contact. The exosphere and bow shock are scaled to the 
same density. The density in the bow shock can be seen to fall off 
rapidly away from the nose of the shock. 

The model parameters reported here are fairly tightly 
constrained. For example, setting 9sh = 40° results in 
too large of a projected area of the bow shock and re¬ 
quires the density to become even more highly concen¬ 
trated at the nose. Similarly, smaller values of a do not 
concentrate the material enough at the nose and the ab¬ 
sorption values are too large. The standoff distance must 
also increase to allow the shock to exit the stellar disk 
at t — tmid = —70 minutes. We note that the values 
of Tm = 12.75 Rp and 9sh = 15° derived in our fa¬ 
vored model are similar to those found by Ben-Jaffel & 
Ballester (2013) (16.7 Rp and 10-30°, respectively) for 
their bow shock model of a pre-transit C II absorption 
signal. However, a value of = 16.7 Rp is too large for 
our data since this would imply egress of the bow shock 
from the stellar disk earlier than observed, assuming a 
similar value of 9sh- Variability related to the stellar 
wind and/or corona, which has likely been observed to 
cause variations in the exospheric absorption signature 
from HD 189733b (Lecavelier des Etangs et al. 2012), 
could account for this difference. Thus we believe our 
values represent a plausible scenario for producing the 
measured absorption. 

The incoming material in a bow shock is not stopped 
by the magnetosphere: momentum is conserved and the 


material flows around the bow (Wilkin 1996). Our model 
does not incorporate the bulk flow of material along the 
bow. The excited hydrogen density we find from our 
model is thus representative of the physical state of the 
gas at the nose of the shock. The rapid falloff of the 
density from the nose implies that the shocked material 
cools as it moves away from the nose and that the in¬ 
crease in temperature and density caused by the shock 
at distances far from the nose is insufficient to produce 
a significant amount of excited hydrogen. 

The large values of and small values of 9eh found 
here and in the study of Ben-Jaffel & Ballester (2013) 
contrast with the relatively small values of (~3.8 Rp) 
and large values of 9sh (~65-70°) predicted by Llama et 
al. (2013) using a full 3D MHD simulation of the stellar 
wind from HD 189733 and assuming a dipolar magnetic 
field strength of 14 G at the pole. The large discrepancies 
suggest that these simulations overestimate the stellar 
wind speed and/or temperature or that the planetary 
magnetic field is underestimated. 

The physical conditions in the bow shock are expected 
to fluctuate as the stellar wind speed, density, and tem¬ 
perature likely vary with time. For example, as the stel¬ 
lar wind speed increases the angle 0^^ will increase, mov¬ 
ing the nose of the shock closer to the star-planet line. 
The low value of 9sh = 15° reported here indicates that 
the planet is moving through a relatively slow moving 
wind, ^30-40 km s“^, or perhaps a hydrostatic corona 
if the true angle is even smaller (Vidotto et al. 2010). 
This may indicate that the planet is moving through a 
quiet or slow region in the stellar wind. Assuming that 
the planetary magnetic field strength remains constant, 
the standoff distance r^. will also change depending on 
the varying thermal and ram pressure from the stellar 
wind. Even if the conditions for a shock to form are not 
met, the incoming stellar wind may still undergo some 
compression and heating as it is slowed by the magneto¬ 
sphere (Lai et al. 2010). Fluctuations in the pre-transit 
signal could be used as a probe of stellar wind variability. 

As a first attempt at using our measurement to esti¬ 
mate properties of HD 189733’s wind, we can use our 
derived model parameters to calculate a mass loss rate 
for the star. If we assume a spherical outflow of con¬ 
stant velocity, the density we estimate at the orbital dis¬ 
tance of HD 189733b translates into a density at the 
stellar surface of po=l-8xl0“^® g cm“^. Assuming a 
wind velocity of 40 km s“^, we estimate a mass loss rate 
of M=2.5xl0^^ g s“^, or ~4xl0“^® M© yr“^. This 
value is ^lOOx lower than that predicted by Llama et 
al. (2013) (4.5x10“^^ Mq yr“^). This suggests that a 
small fraction of the stellar wind material ends up as 
excited hydrogen in the bow shock. 

This estimate is a lower limit: the density from our 
model is the density of the excited hydrogen in the bow 
shock and does not include ground state hydrogen or 
other elements. This mass loss estimate has other obvi¬ 
ous limitations. For example, the stellar wind is likely 
not spherically symmetric. In addition, this estimate is 
essentially a snapshot in time. HD 189733 is an active 
star and likely has a variable wind. Our estimate is also 
lower than the mass loss rate suggested by the X-ray 
flux versus M relationship given in Wood et al. (2014). 
Using an average X-ray flux of Fjc =6.3x10^ taken from 







the NEXXUS2 catalog^, we calculate a mass loss rate 
of M=l.lxlO“^^ Mq yr“^. Although the Wood et al. 
(2014) relationship is fairly uncertain, this value is more 
in line with the Llama et al. (2013) prediction than the 
value estimated from our model. Thus the low mass loss 
rate estimated here should be considered a rough lower 
limit for HD 189733. Future observations will help con¬ 
strain the variability and potentially measure changes in 
the stellar mass loss rate. 

6 .2. Exosphere model 

The exosphere model consists of an isothermal gas in 
hydrostatic equilibrium above the planet. We calculate 
the density in the exosphere out to lOiLe^o where iLexo is 
the exospheric scale height given by H^xo = kbT/ (rrihg) 
where is Boltzmann’s constant, T is the exospheric 
temperature, and g is the acceleration due to gravity at 
the planetary surface. The exobase is taken to be at 
i?p=1.0. In order to roughly match the in-transit ab¬ 
sorption, we find Hexo = 0.04 Rp and T^xo = 7,000 K 
with a base number density of ng = 300 cm“^. This sim¬ 
ple model produces systematically low estimates for the 
H/? absorption. This likely indicates that the exosphere 
does not follow a simple hydrostatic pressure law and 
a more realistic treatment is necessary to more closely 
match the line ratios. We note that the exospheric pa¬ 
rameters derived here and very similar to those derived 
by Christie et al. (2013) who give a detailed treatment 
of the n = 2 hydrogen level population in the exosphere 
of HD 189733b. 

7. ESTIMATING THE PLANETARY MAGNETIC FIELD 
STRENGTH 

Our measurements can be used to estimate the plan¬ 
etary magnetic field strength by considering some limit¬ 
ing cases of pressure balance between the interplanetary 
material (i.e., stellar wind or corona) and the planetary 
magnetosphere. Rewriting equation (9) from Llama et 
al. (2013) we see that the magnetic field strength of the 
planet at the pole is 

Bp = 2 (^^^ [8^(p^Au^ + P„)+B2]1/2 (4) 

where is the mass density of the stellar wind or corona, 
Auw is the relative velocity of the planet and the wind, 
Pw is the thermal pressure of the wind, and is the in¬ 
terplanetary magnetic field strength (i.e., the stellar mag¬ 
netic field strength at the orbital radius of the planet). 
All of the pressure terms are the values at the orbital ra¬ 
dius of the planet. Equation 4 neglects any pressure from 
the planetary atmosphere or wind. The equatorial mag¬ 
netic field strength, assuming the field is purely dipolar, 
is one-half the polar value. 

We adopt our favored model parameters and assume 
that the bow shock is mediated by a magnetosphere. If 
we neglect the pressure due to the stellar magnetic field, 
the dipolar magnetic field strength of the planet can be 
estimated by assuming pressure balance between the stel¬ 
lar wind and the planetary magnetosphere (Vidotto et al. 

^ http://www.hs.uni-hamburg.de/DE/For/Gal/Xgroup/nexxus/ 


2010; Llama et al. 2013). By neglecting the stellar mag¬ 
netic pressure we are essentially estimating a lower limit 
for the field since any addition to the pressure will require 
a larger planetary field strength for the inferred stand¬ 
off distance. We choose a stellar wind temperature of 
Tyj = 1 X 10® K which is roughly consistent with a wind 
of velocity of <100 km s“^ at 0.03 AU (Matsakos et al. 
2015). The wind is assumed to be traveling radially at 40 
km s“^. The azimuthal wind velocity is neglected. The 
stellar wind is taken to be an ideal gas composed purely 
of protons to provide an estimate for . We also assume 
that the material in the shock has been compressed by 
a factor of 4, i.e., the shock is adiabatic, resulting in a 
stellar wind density of /9u,=2.5x 10“^® g cm“®. We note 
that this density is very similar to the stellar wind den¬ 
sity derived by Bourrier & Lecavelier des Etangs (2013) 
(^3x10“^® g cm“^) to reproduce observed La absorp¬ 
tion in HD 189733b. 

Given these assumptions we estimate an equatorial 
magnetic field strength of B^q = 28 G for HD 189733b, 
similar to the upper limit of 24 G derived by Vidotto et 
al. (2010) for WASP12-b. The estimated field strength 
is ~7 times larger than Jupiter’s equatorial field (~4 G) 
and is driven by the large value of r^n- The scaling rela¬ 
tion of Reiners et al. (2009) is used by Reiners & Chris¬ 
tensen (2010) to estimate the magnetic field strength of 
HD 189733b at ^14 G. The factor of two difference be¬ 
tween our estimates and those from the scaling relation 
indicate a factor of ~8 difference in the internal heat flux 
used in the scaling relation estimate. The additional heat 
source suggested by our estimates could have significant 
implications for models of hot Jupiter interiors. 

While we have focused on a planetary bow shock as the 
mechanism for producing the pre-transit absorption, al¬ 
ternate scenarios exist for producing similar signals. The 
effects of charge exchange between a hot stellar wind and 
a planetary wind are investigated in Tremblin & Chiang 
(2013). Their model, which has a similar geometry to the 
bow shock approximation, produces a layer of hot, dense 
neutral hydrogen at a similar standoff distance (^10 Rp) 
as found in our model and shows decent agreement with 
Lyman a absorption measurements. Their model likely 
predicts too large of an excited hydrogen density over too 
large a volume to match the low absorption values found 
here. However, a more detailed investigation of excited 
hydrogen absorption in this context would be informa¬ 
tive. Indeed, charge exchange may be occurring in the 
bow shock between stellar wind protons and neutral hy¬ 
drogen that has escaped from the planet, contributing to 
the n=2 hydrogen density. Transiting accretion streams 
from the planet, which is overflowing its Roche lobe, to 
the star can produce sufficient column density and cover¬ 
age of the stellar disk to result in pre-transit absorption 
(Lai et al. 2010). This scenario, however, should produce 
red-shifted absorption which is not observed in our data. 
In fact, we find no absorption centroids shifted by more 
than ±5 km s“^ from line center. This suggests that the 
excited hydrogen has a symmetric bulk velocity along the 
line of sight. 

8 . SUMMARY AND CONCLUSIONS 

We have presented our detection of resolved pre-transit 
line absorption around HD 189733b in the first three 
members of the hydrogen Balmer series. A bow shock 
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model is able to reproduce the important features of the 
absorption time series. The physical parameters derived 
from our model suggest an equatorial planetary magnetic 
field strength of Beq = 28 G. This large value is driven 
mainly by the large standoff distance of the magneto¬ 
sphere rm = 12.75 Rp, a parameter that is well con¬ 
strained by the model. For the derived standoff distance, 
this field strength is likely a lower limit due to the mul¬ 
tiple simplifying assumptions employed in the estimate. 
More detailed models which include the stellar magnetic 
field pressure and a full treatment of the stellar wind 
will prove insightful. The expected variability of the 
bow shock signal can also be tested by observing mul¬ 
tiple transits across many epochs. Observations of this 
type could be useful probes of stellar wind variability for 
HJ hosts. 

The magnetic field of a HJ plays an important role in 
determining the planetary mass loss rate: the stronger 
the magnetic field, the lower the outflow rate (Owen & 
Adams 2014). Planetary magnetic fields also provide 
a potential source of energy for SPIs via reconnection 
events with the stellar corona (Cuntz et al. 2000). This 
energy flux scales as Bpi (Cuntz et al. 2000; Lanza 2013) 
so larger magnetic field strengths can potentially induce 
more energetic SPIs. If the field strengths estimated 
above are common for HJ exoplanets, this regime needs 
to be more fully explored by mass loss and SPI models 
which have generally considered field strengths similar to 
or less than that of Jupiter (Lanza 2013; Trammell et al. 
2011; Strugarek et al. 2014). 

Planetary magnetospheres are a revealing diagnostic 
of planetary interiors and the first line of defense against 
energetic stellar wind particles which can have a dev¬ 
astating impact on planetary atmospheres and life near 
the surface of the planet. While considered to be a crit¬ 
ical property in understanding a planet, its evolution, 
and habitability, the prospects of measuring an exoplan¬ 
etary magnetic field were thought to be remote. How¬ 
ever, if independent measurements can be made of the 
planet hosts’ stellar wind, observations such as those pre¬ 
sented here may provide important constraints on exo¬ 
planetary magnetic fields. The large field strength we 
estimate implies that lower magnetic field strengths, as 
found in our solar system, may not be ubiquitous among 
other exoplanet systems. A population study of magnetic 
field properties of exoplanets that display bow shock-like 
structures will complement measurements of solar sys¬ 
tem magnetospheres and be used to better understand 
planetary interiors and the intertwined relationship of 
magnetic fields, planetary atmospheres, and habitability. 
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